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From QCD to heavy ion collisions 
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This talk will discuss how heavy ion experiments, when moving from SPS (10 + 10 GeV) 
to RHIC (100+100 GeV) and to LHC (2750+2750 GeV), will enter a new domain of QCD 
in which the production of even large px gluons is so abundant that it is simultaneously 
perturbative but such that the phase space density of gluons is saturated. The saturation 
scale p sat is estimated, quantitative numbers for the initial production of gluons at the 
t^j- ! LHC are given and options for their subsequent evolution are outlined. For parametrically 
| large nuclei and energies, classical field methods will be applicable. 

O 

ON • 1. Perturbative calculability in heavy ion collisions 
0\ 

The first collisions at RHIC are expected to take place just during this PANIC confer- 
ee ence and the outcome for, say, average charged multiplicity in a central collision will be 
Ph! known by the time the Proceedings are out. However, in spite of the fact that QCD is 
a perfectly well known theory we cannot with reliability predict this outcome - let alone 
the properties of QCD plasma produced in these collisions. 



In perturbative QCD reliable predictions can be made for infrared safe quantities when 
there is a hard scale Q in the problem, A prototype is inclusive large jet production, 
in which the hard (the parton-parton collision) and the soft (the structure functions at 
scale Q) parts of the problem can be factorised: 
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This schematic diagram is explicitly drawn for a nucleus-nucleus collision, though calori- 
metric jet measurements will need very large Et due to the immense hadronic Et back- 
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ground. Note also that the unobserved X's will in an average event contain a huge number 
of further 2^2 collisions. 

Consider now the situation as a function of px (the scale is for LHC): 

12 100 GeV 

i 1 1 1 ^ pT 

Classical field <^ Perturbation theory, &(pt) small, 

equations ^ minijets error ~ 10% 

Psat 

When one reduces the magnitude of pt, the cross section grows and the error likewise. 
In this perturbative minijet picture one can compute the number of gluons with p? > p 
and the Et they carry [Q-Q. We can continue until po reaches the saturation limit p sat , 
to be defined below. Then the gluonic subsystem becomes very dense and new physics 
enters For px <C p psa ,t the phase space occupation number becomes 3> 1. Then 

this new physics may be effectively described by classical Yang-Mills equations 0, the 
initial conditions of which are given by an ensemble in colour space [§J-0- 

It is now crucial that for large nuclei and large energies the gluonic subsystem may 
become so dense that p sat is in the perturbative region, >2 GeV. 

Numbers for p T 3> p sa t can easily be produced using standard perturbative techniques. 
Classical field ideas are not yet sufficiently developed to permit the same for Aq CD <C 
Pt "C Psat and may ultimately require parametrically large nuclei (A 1 / 3 ^> 1) and energies. 
To have some quidance from theory let us use the assumption already advocated in [[J: 
estimate the average event by taking pt = p sa t- We are mainly interested in the local 
energy density or Et (per unit rapidity): gluons with pt 3> p s &t carry lots of Et but are 
very rare, whereas gluons with pt <C p sa ,t are numerous but carry little Et- 

2. Minijet approach to the creation of early little bang 

During the pre-little bang era we have two Lorentz contracted nuclei approaching each 
other. The system has zero temperature and zero entropy. The Lorentz contracted 
diameters, applicable to the valence quark parts of the wave functions, are 

^ = lfm SPS, (1) 

7 

= 0.1 fm RHIC, 
= 0.005 fm LHC. 

Already this naive fact emphasises the great qualitative difference between the energy 
domains. The little bang universe now is created by the nuclear collision, treated in 
the standard picture as a collision of two clouds of gluons (which dominate over quarks 
and antiquarks) with longitudinal momentum distributions ~ Axg(x,p T ). The entropy 
content of this universe can be bounded from below by computing the number of gluons 
with pt > Po = 2 GeV. The dominant longitudinal momentum fractions are (when gluons 
from the two nuclei have equal pi = Pt) 

x ~ — =0.2 SPS, 
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0.02 RHIC, 



0.001 LHC. 



This is a second naive way of seeing the importance of increasing energy (keeping p 
fixed): LHC can fully profit from the enhancement of gluons observed at small x [ |T0|j . 

The computation of the number Naa of gluons produced in a central A+A collision 
is quite straightworward ||. What is important for a qualitative discussion is that this 
number scales ~ v4 4//3 . A hard cross section scales ~ A 2 , but since one wants the number 
per inelastic collision, this has to be divided by ttR\ ~ A 2 ^: 



N, 



AA 



A 2 a(p T > Po) _ .4/3 J_ 



71 R 2 



1 ' 

Po 



(2) 



where a is a 2 — » 2 cross section. Computed values per unit rapidity at y — are as 
follows (shadowing is included, K = 1.5): 

<N> of minijets in central rapidity unit 
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The result is reliable as long as po is large enough so that perturbation theory is valid 
and as long as the independent collision picture can be assumed to be valid. The latter 
becomes questionable when p < p sat , the saturation limit. To define p sat , consider the 
transverse projection of the nucleus: 
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One can estimate the each gluon occupies the transverse area ir/p? so that the gluon 
density saturates if 

Naa x n/pl t > ttR 2 a (3) 

from which it follows that p sat ~ A 1 ^. This saturation limit is also shown in Fig.?? and 
one sees that 

p sat = 0.65 GeV SPS, (4) 
= 1.1 GeV RHIC, 
= 2.0 GeV LHC. 

Extending the computation to include partons with p T > p sat one has, in one unit of 
rapidity: 

SPS: 600 gluons with p T > 0.65 GeV, 

RHIC: 1100 gluons + 110 q + 80 q with p T > 1.1 GeV, 

LHC: 4300 gluons + 200 q + 190 q with p T > 2.0 GeV. 
The corresponding .E^- values are 

E T = 400 GeV SPS, (5) 
= 2500 GeV RHIC, 
= 12000 GeV LHC. 

Remember that these are the initial values at r = Tj = l/p sa t- 

It is, of course, questionable to extend the computation to the 1 GeV region, but the 
numbers are not unreasonable. However, due to the (relatively) small Lorentz contraction 
at SPS it is not clear how to convert to coordinate space there. At LHC the clocks can 
be started very accurately and one can say that at the LHC there will be at the time 
Tj = l/p sa t = 0.1 fm after the collision about 4000 gluons with p T > p sat in one unit of 
rapidity: 



Geometrically, the volume of the little bang universe at the time of its creation is 
V = nR A -k n Ay = 7iR A /p sat = 13 fm 3 for A = Pb. This implies that the initial energy 
and number densities of gluons at the LHC are 

€i = 970 ftn^' n * = ^ (6) 
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3. Hubble expansion of the early little bang 



Now comes the crucial question: what then? We have in usual QCD perturbation 
theory computed the composition of the > Psat = 2 GeV component of the system 
at the time l/p sa t = 0.1 fm. Invoking the assumption of ||, this is the initial state of 
an average event. But this is not yet something comparable with experiment. How does 
this subset behave for r > l/p sa t until conversion to hadrons? Of course, nobody knows 
definitely until a systematic non-perturbative computation in QCD has been carried out. 
Thus one has to make assumptions and these will be tested by experiment. 

In this case it is quite natural to assume that the little bang universe is thermalised at 
creation at r = 0.1 fm (note that chemically (anti)quarks are not in equilibrium). Then 
it expands conserving the total entropy S = 3.60iV w 15000. If we go to the center, the 
expansion is purely longitudinal, S = sV ~ st ~ T 3 r = constant, s ~ n ~ 1/r and 
T ~ 1/r 1 / 3 (compare T ~ l/t 1 / 2 (~ 1/i 2 / 3 ) for the usual radiation (matter) dominated 
universe). However, due to pdV work, energy density decreases faster |Tl| : 
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energy density vs. proper time 
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Thus the initial and final entropies are the same. The final particles are dominantly 



pions, for which also S 

1/3 

E- 
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iVpionJ Et behaves quite differently: 



(8) 



Thus the Et goes down from 12000 GeV to about 2000 GeV. This plasma really does 
work! Stated in another way: the energy per particle went down from about 3 GeV to 
about 0.5 GeV. These predictions for Pb+Pb at LHC can easily be extended to any A 
and s. 

Any dissipative effects will increase both the multiplicity and Et- The extreme is 
free streaming: Et constant. In this limit the multiplicity would be gigantic, about 
-Er/O.SGeV. If one could predict the initial parameters with a 10% accuracy, the mea- 
surement of final quantities would be a measurement of the degree of thermalisation. Or, 
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if one knew that the expansion is adiabatic, they would be a measurement of the initial 



conditions [O. 



4. Classical fields 

In the minijet computation above we basically used the 2^2 diagram (x denotes a 
large px external line): 




but what about bremsstrahlung diagrams, like 





This diagram predicts |H| a certain number of gluons around y = 0, but how to extend 



this computation to nuclear collisions? This is not as straightforward as in the minijet 
case: the px of the perturbative gluon is balanced by a "beam jet". The new ideas 
developed [§[0, |II|-[|I7J amount to computing classical radiation of gluons by two colour 



currents, formed by the valence quarks and the fastest gluons of the colliding nuclei, on 
the forward and backward light cones. 

The motivation for the use of the classical field approximation is precisely that in 
the saturation limit the occupation number of any of the field modes is large and the 
classical field approximation should be appropriate. The situation is analoguous to that 
of computing the rate ( "the sphaleron rate" ) of the baryon number violating reactions in 
electroweak matter ||18|| . These reactions also involve nonperturbative large field modes 
and their rate can only be computed numerically. 

The essential parts of the classical field computation are the formulation of the equations 
of motion ("effective theory") and of the initial conditions. In the electroweak case, 
the standard model at finite T, the theory is weakly coupled and there is a perfectly 
controllable perturbative way |19| of computing the effective theory, which is a three- 
dimensional purely Euclidian SU(2)xU(l) gauge + fundamental (doublet) Higgs theory. 
The initial conditions are set by thermalising the system at some temperature, which 
is numerically straightforward. After that the classical field equations are integrated 
numerically in real time starting from some thermal field configuration and the sphaleron 
rate is computed as a real time correlation function of the topological susceptibility. 
Hereby one has learnt that the dynamics of the infrared small-momentum modes one 
is interested in is essentially affected by their interactions with the hard large-momentum 
(~ T) modes. When this is correctly taken into account, the calculation of the sphaleron 
rate, which has no perturbative contribution, should be under control [[20]. This is not 



yet the case with quantities like viscosities, which also have a perturbative contribution. 
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The situation is much more complicated in the application of the classical field approach 
to QCD in heavy ion collisions. Here there is no "top down" systematic way of deriving 
the effective theory, even in a theoretical weakly coupled limit and one has to proceed 
phenomenologically ||. Let us focus our attention at the slice of matter at rest in the 
CMS. Then our space is again three-dimensional, but in this case Minkowskian with 
time and two spatial coordinates. The effective theory now is an SU(3) gauge + adjoint 
(octet) Higgs theory |17]], but there is no known systematic way to derive its coefficients. 
To see why the theory is of this form, take the gluon part of the QCD action and include 



only longitudinal boost invariant configurations, independent of p 
implies an infinite energy. The action then becomes 



drd x 



T 1 

4 2t K 



u S 



) a {D^ 



tanh (z/t), which 



(9) 



where p = 0, 1, 2, A±, are the transverse components of the gluon field, a is the 
component of the gluon field and the gauge has been chosen as A a T = 0. This part 
of the action is singular at r = 0, on the light cone. This singularity is regulated by 
the dynamics of the nuclear interaction at r = 0, by the initial condition. One collision 
is represented by one random charge distribution p^xi, £2) of the two colliding nuclei 
(i = 1,2), where each of the two p's is drawn from a, say, Gaussian distribution of some 
width p. The p's give the initial conditions for the fields ^(0, x±, X2), a (O, x±, 22) and 
the equations of motion give their further evolution in r. Various physical quantities, 
like the number of bremsstrahlung gluons, can then be computed by averaging over many 
different evolutions corresponding to different initial conditions drawn from the Gaussian 
ensemble. A first attempt to carry out this program (for SU(2)) is in [|T7j . 



The crucial parameter is the average transverse colour charge density p 2 , which for the 
valence quarks of a nucleus is given by 



(iV 2 -l) A 



2iV 2 



(10) 



but to which also the (dominating) contribution of the fast gluons has to be added |L5 
The computation then is formally valid for 



Aqcd < Pt < P, 



where the lower limit expresses the fact that one is in the weakly coupled domain (a s is 
parametrically small) and the upper that one simultaneously is in the saturation domain, 
occupation numbers are large. Parametrically, like p sat , p scales ~ A 1 / 6 . However, numer- 
ically at LHC energies, p ~ 0.8 GeV |I5| , while p sat was about 2.0 GeV. Converted to time 
units this means that at LHC the bremsstrahlung quanta thus will materialise at times 
0.25 fm < r < 1 fm, while the perturbative quanta materialised at r pa O.lfm. For practical 
applications, the window of applicability of the classical field approach thus seems very 
small, but conceptually it deals with an entirely new domain of QCD and thus studies 
with parametrically large A, 1/g and \/s are very important. 

However, even if one could solve the problem for all the scales px ^> Aqcd, implying 
the computation of the initial creation of the system at all times -C 1 fm, there still 
remains the problem of all smaller p^-scales, the entire further evolution of the system, 
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its expansion, conversion to the hadronic phase and ultimate decoupling. Here one in 
any case will have to resort to phenomenological ideas, which later will be checked by 
experiments. 



5. Conclusions 

Heavy ion experiments at RHIC and LHC will enter a new domain of QCD, in which 
the production of even large px gluons is so abundant that it can be simultaneously 
perturbative but such that that the phase space density of gluons is saturated. This talk 
has discussed precisely what numbers of gluons can be expected under these conditions to 
be initially (at 0.1 fm) produced at the LHC and what options there are for the subsequent 
behaviour of this system. Whatever the detailed validity of these numbers is, one can in 
any case with confidence look forward to abundant production of QCD plasma at RHIC 
and especially at the LHC 
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